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Abstract. We present model calculations for the expected surface density of Ly-a emitting primeval galaxies 
(PGs) at high redshifts. We assume that elliptical galaxies and bulges of spiral galaxies (= spheroids) formed 
early in the universe and that the Ly-a emitting PGs are these spheroids during their first burst of star formation 
at high redshift. One of the main assumptions of the models is that the Ly-a bright phase of this first starburst in 
the spheroids is confined to a short period after its onset due to rapid formation of dust. The models do not only 
explain the failure of early surveys for Ly-a emitting PGs but are also consistent with the limits of new surveys 
(e.g. the Calar Alto Deep Imaging Survey - CADIS). At faint detection limits Sj: m < 10~ 20 W/m 2 the surface 
density of Ly-a emitters is expected to vary only weakly in the redshift range between z — 3 and z = 6 with values 
> 10 3 /deg 2 /Az = 0.1 reaching its maximum at zo w 4. At shallower detection limits, S-\: m > 3 x 10 _20 W/m 2 
the surface density of high-z Ly-a emitters is expected to be a steep function of redshift and detection limit. 
This explains the low success in finding bright Ly-a galaxies at z > 5. We demonstrate how the observed surface 
densities of Ly-a emitting PGs derived from recent surveys constrain the parameters of our models. Finally, we 
discuss the possibility that two Ly-a bright phases occur in the formation process of galaxies: An initial - primeval 
- phase in which dust is virtually non-existant, and a later secondary phase in which strong galactic winds as 
observed in some Lyman break galaxies facilitate the escape of Ly-a photons after dust has already been formed. 
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1. Introduction 

The detection of the ancestors of large present day galaxies 
(like our Milky Way) during their first phase of violent star 
formation is currently one of the great challenges for obser- 
vational cosmology. Commonly, these objects are referred 
to as primeval galaxies (PGs). Finding them in substantial 
numbers and over a sufficiently broad range of luminosi- 
ties would provide us with direct insight into the epoch 
of galaxy formation in the young universe. Since current 
wisdom places the PG phase of our Galaxy between red- 
shifts z = 5 and z = 10 one should ultimately aim to 
establish the luminosity function of primeval galaxies and 
their evolution at several redshifts within this range. 

It was not long ago that such observational program 
would have seemed audacious. But with half a dozen tele- 
scopes of the 8-10m class in operation, the detection of star 
formation rates of ten Mq / year appears feasible even at 
z > 5. Indeed, the last few years have witnessed an enor- 
mous progress in identifying young galaxies at very high 
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redshifts (see Dey et aL [13551 Weymann etal. IT3351 Hu 
etal. [13551 IT3331I2UU21 I2TJU51 Rhoads etal. I2oTJ5|i . 

Various techniques have been used to search for 
galaxies at high redshifts. Currently the most successful 
metho d, intro duced b y Steidel et al. <(T3521 113551 11996a! 
11996b! 11998a! I1998bjl uses the Lyman break and the 
flat spectral energy distribution blue-wards of the Lyman 
break as a signature of very distant, young star- forming 
galaxies. Hundreds of galaxies found in this way have 
been spectroscopically confirmed to be young star forming 
galaxies at redshifts z ~ 3 (Steidel etal. 1199 6a). In addi- 
tion several dozens of them could be detected at z ~ 4 
(Steidel et al. I199"5jl . The number density and clustering 
properties of these Lyman break galaxies are consistent 
with them being the central galaxies of the most mas- 
sive dark matter halos present at z ~ 3 (Mo etal. 1999). 
Although the Lyman break galaxies seem to be relatively 
young star forming galaxies, they cannot represent the 
population of primeval galaxies in the sense defined above: 
The strong metal absorption lines in their spectra indi- 
cate that their current star formation period must have 
been preceded by an earlier star burst. The massive stars 
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formed have already substantially enriched the interstellar 
medium in these systems. Moreover, both the UV contin- 
uum slope and the Balmer lines (Pettini et aZ l!998|l suggest 
that the Lyman-break galaxies already contain significant 
amounts of dust. Therefore, it is now widely accepted that 
the star forming rates inferred from the UV continua of 
Lyman-break galaxies have to be corrected upwards by 
a factor of 2 to 5. The presence of dust also explains 
why only about 1/3 of the Lyman-break galaxies exhibit a 
strong Ly-a line: In typical star forming regions every Ly- 
a photon emitted by the hot OB stars will undergo dozens 
of multiple resonant scatterings before leaving the region 
towards the observer. Thus even small amounts of dust 
could quench the escaping Ly-a emission considerably. 

In contrast to the Lyman-break galaxies, all galaxies 
known at z > 5 exhibit a very strong Ly-a emission line 
(rest frame equivalent widths > 5 nm) . This is exactly the 
spectral signature we expect for the first few hundred mil- 
lion years after the onset of a violent burst of massive star 
formation and before the newly produced metals could be 
recycled into the cold phase of the inter-stellar medium. 
Note that at z = 3 there still exists a large population 
of Ly-a bright galaxies (Shapley etal. 120031 Kudritzki 
et al. I2000fl but their low star forming rates indicate that 
they represent a population of smaller galaxies in which 
the "trigger density" has been reached later than in the 
Lyman-break systems. The weak continua of Ly-a bright 
galaxies make them hard to find by color-break techniques, 
but the strong emission lines should easily be detected in 
narrow-band searches for emission line objects. 

Although early attempts to detect Ly-a bright galax- 
ies at z 3 failed (Pritchet 1X9948) , we have witnessed a 
breakthrough in finding these objects at redshifts between 
3 and 5 as sufficiently deep detection limits (line fluxes of 
a few x 10 -20 Wm" ! ) have been reached routinely in the 
last years (Hu et al. ll^U^^OHl Rhoads et al. I3UU31 
Maier etal. .2003). The most distant known object in the 
universe has been found in this way (Hu et al. 2002 ) . But 
still, the number of Ly-a bright galaxies at z > 5 which 
have been found in systematic surveys is very limited, and 
larger samples are required to draw any firm conclusions 
about the epoch of galaxy formation. 

For both the interpretation of the results of present 
narrow-band searches for Ly-a bright galaxies and the op- 
timum design of future surveys, it is essential to estimate 
the expected abundance of these objects under reasonable 
assumptions about the cosmological parameters and the 
history of galaxy formation. 

Here we present a phenomenological model to pre- 
dict the expected surface density of Ly-a bright PGs 
at high redshifts based on ideas which have already 
been sketched in Thommes & Meisenheimer l|1995|l and 
Thommes fiWty . 

In principle, the prediction of the abundance of Ly-a 
bright galaxies at high redshifts can be obtained in two 
ways: One way starts from a primeval density field in 
the early universe and follows the collapse of dark matter 
haloes by N-body simulations or with the Press Schechter 



formalism. Adding in baryonic matter and a reasonable 
star formation scenario in combination with a treatment 
of dust formation and distribution could then predict the 
abundance of star forming haloes and their star formation 
rate which - under the assumption of an initial mass func- 
tion (IMF) - could be converted into a prediction of the 
number density of Ly-a bright galaxies above a certain 
detection limit. Haiman and Spaans (1999) present calcu- 
lations along this line. However, it is not trivial to scale 
such an ah initio approach to the observed abundance of 
galaxies in the local universe. 

Therefore, we pursue the second way in which one tries 
to extrapolate the local luminosity function of galaxies and 
their stellar content back into the past. This way has been 
pioneered by Meier (1976) and further explored by Baron 
& White l|1987[l leading to predicted surface densities of 
Ly-a bright PGs at z ~ 5 between 1.4 PGs/D'/(Az = 1) 
(for q = 0.5) and 0.05 PGs/D'/(Az = 1) (q = 0.0) for a 
survey limit Su m > 10 -19 Wm -2 . Such a high abundance 
has been clearly ruled out by the narrow-band surveys 
carried out by Thompson et al. ( 1995a) and more recently 
by Hu etal. (1998). This has triggered several new efforts 
to search wider fields to deeper limits (e.g. the Calar Alto 
Deep Imaging Survey - CADIS see Meisenheimer etal. 
irMTllTMSl the Large Angle Lyman Alpha survey - LALA, 
see Rhoads et al. 120001 . We have identified two main points 
which Baron & White (Qt)87) did not take into account: (1) 
The Ly-a bright phase of PGs might be rather short due to 
rapid dust formation, and (2) galaxies show a substantial 
age spread and therefore did not form or start simultane- 
ously with their first star formation. As we will see, both 
points tend to reduce the expected number densities so 
that the apparent contradiction between observations and 
predictions disappears. 

In the present paper we will describe the basic assump- 
tions of our models and discuss how choices of the cosmo- 
logical parameters and the history of galaxy formation and 
of primeval star formation would influence the observable 
number density for a very broad range of search redshifts 
3 < z < 13 and detection limits between 10 -19 and 10 -21 
Wm~ 2 . In a forthcoming paper (Meisenheimer, Thommes 
and Maier 2004, in the following referred to as paper II) we 
will try to combine all available survey results to constrain 
the free model parameters even further and thus provide 
a much more confined set of predictions which could be 
used as bench mark for future surveys. 

The present paper is structured in the following way: 
In section 2 we briefly describe the narrow-band imaging 
technique for detecting Ly-a bright galaxies. In section 
3 we describe the principle assumptions, parameters and 
functions of our model. In section 4 we give a very simpli- 
fied and transparent version of our model and demonstrate 
why the earlier predictions by Baron & White (1987) were 
much too optimistic. Section 5 explores the full range of 
model parameters in order to identify those which will 
most critically affect the predicted number density of Ly- 
a bright primeval galaxies. In section 6, we summarize 
the generic results of our model. Section 7 discusses how 
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our results are useful in designing optimum surveys and 
gives a first account of how well the model agrees with 
the results of present surveys. This issue will be detailed 
in the subsequent paper II, in which we will adjust the 
free model parameters as close as possible to the results 
of all available emission line surveys. This will constrain 
the range of free parameters even further. 

2. Narrow band imaging technique to search for 
high-redshift Ly-a galaxies 

Any emission line survey must aim to map the three- 
dimensional phase space of objects (a, 8, \ bs), where 
a, 5 are the positions on the sky and A fc s = (1 + z)\ e i 
is the observed wavelength of the emission line, onto the 
two-dimensional detector in an optimum way (A e ; = rest- 
frame wavelength of the emission line). Standard obser- 
vational techniques for emission line surveys are reviewed 
by Pritchet lfT§M|l . 

We concentrate on the narrow-band imaging tech- 
nique (see e.g. Hippelein et al. 1271(131 Maier etal. 2003 
Meisenheimer etal. 1997): Here a rather narrow wave- 
length range AA is selected by using a narrow-band filter 
or an imaging Fabry-Perot-Interferometer, while (a, S) is 
directly mapped onto the detector coordinates (x, y). This 
technique provides two main advantages for the search of 
Ly-a galaxies at z > 5: (a) Since the entire field of view 
Auj is mapped onto the detector, two-dimensional meth- 
ods can be used for background determination and source 
extraction. In general, they perform much better than the 
one-dimensional methods used in analyzing long-slit spec- 
tra, (b) One can place SX n such that it falls into wave- 
length regimes of low and smooth night sky emission. This 
is of high importance when searching for Ly-a emission at 
z > 4.9 where about 75% of the wavelength range is made 
rather useless by strong OH-lines in the night sky. 

For the following model predictions of the abundance 
of Ly-a galaxies we assume Az = AA/Alj, q = 0.1 
(Aij /Q =121.57 nm) which is typical for the narrow-band 
technique. Three redshift intervals are selected such that 
the Ly-a line falls into the best night sky windows around 
A ~ 705, 815, and 920 nm (that is z ~ 4.8, 5.7, and 6.6). In 
addition, we present predictions for z = 3.5 (A = 530 nm, 
V band), z = 9.3 (A = 1250nm, J band), and z = 12.6 
(A = 1650 nm, H band). 

3. The basic formalism 

We want to calculate the surface number density of Ly- 
a emitting PGs on the sky per solid angle Auj which 
have detectable Ly-a fluxes greater than a certain flux 
limit Su m and which have redshifts zq in an interval 
[zq — Az/2, zq + Az/2]. In a narrow band search for Ly-a 
emitting objects, zq is given by the central wavelength Ao 
of the narrow band filter and Az by the band width of the 
filter AA (Az = A\/\ Lya ; A iya =121.57 nm). Here Auj 
refers to the area on the sky which is covered by the entire 
survey. Aw, Az and zq define a certain comoving volume 



AV C of the universe. The intrinsic Ly-a flux of PGs at zq 
which produce observable fluxes of Su m is given by 

L min (zo) = 47r Sum d L (z ) 2 (1) 

with the luminosity distance (Il(zo)- We have to calcu- 
late the number of PGs in the volume AV C with Ll voi > 
L m in{zo) at the epoch to = t(zo). This number depends 
on: 

(i) The galaxy formation history, which determines the 
number of galaxies in the volume AV C which are just 
in the process of forming and therefore might be in 
their Ly-a bright phase at the redshift Zq. 

(ii) The evolution of the Ly-a luminosity of the PGs as a 
function of time and mass of the PG. 

We assume that the stars in ellipticals and bulges (both 
called spheroids) precede the formation of stars in disks 
which form out of gas which accretes around the spheroids. 
Our hypothesis is that most Ly-a emitting PGs will be 
these spheroids during their first burst of star formation 
at high redshift (but also refer to our discussion in section 
7.4). We assume that such a "proto-spheroid" starts to 
shine in Ly-a as soon as the first star formation sets in 
and that the Ly-a luminosity is proportional to the star 
formation rate (SFR) as long as dust absorption is neg- 
ligible. Furthermore, it is reasonable to assume that the 
SFR of the first star burst is proportional to the bary- 
onic mass of the PG. Thus, the Ly-a luminosity of a PG 
should be proportional to its baryonic mass which we as- 
sume to be proportional to the total mass of the dark 
matter halo surrounding the spheroid. If the SFR is not 
constant but changes with time, so too does the Ly-a lu- 
minosity. However, after some time other effects like dust 
absorption may quench the Ly-a luminosity of the PGs, so 
the proportionality to the SFR is no longer valid. We de- 
scribe this time dependence of the Ly-a luminosity with 
a function f(t — t s ) where t s is the epoch at which the 
galaxy starts to shine in Ly-a (see Fig. Q] we call t s the 
ignition time): 

L ( Ly b a \t-t s )=k Lya M b f(t-t s ) 

(2) 

with f( t -Q l=° n 

Jy ' \ > if t > t s 

khya gives the factor of proportionality between SFR and 
Ly-a flux in a dust-free medium. Mb denotes the baryonic 
mass of the PG which we assume to be proportional to 
the total mass M of the PG. L^ b J(t — t s ) gives the Ly-a 
luminosity of a galaxy with mass Mb at the epoch t (time 
since the big bang), which started to shine in Ly-a at the 
epoch t s . Because of the assumed proportionality between 
the Ly-a flux and the PG mass M, there exists a lower 
limit M m i n for the PGs observed at redshift zq, which can 
get luminous enough in Ly-a that they have an observable 
flux > Sii m . PGs with M < M m i n remain always too faint 
in Ly-a to be detectable above the survey flux limit Su m . 
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Fig. 1. The constraint t s has to fulfill so that the Ly-a 
flux is detectable at zq. The big black point marks the 
minimum flux detectable at zq — z(to) with the detec- 
tion limit Sii m . In order to be detectable in Ly-a, the 
ignition epoch t s for the Ly-a emission has to be in the 
interval [t Q — Ai 2 ,i — Ati]. This interval is marked by 
the thick vertical line. The short dashed, solid and long 
dashed curves show the possible Ly-a luminosity of PGs 
as a function of time with different ignition times t s . The 
short dashed curve shows the time evolution of the Ly- 
a luminosity of a PG with an t s at the lower bound of 
the relevant time interval. At to the Ly-a luminosity of 
such a PG is already very week e.g. because of ongoing 
dust formation. The solid curve corresponds to a PG with 
an ignition time t s centered in the relevant time interval. 
It reaches its maximal Ly-a luminosity at the epoch to 
corresponding to the observed redshift zq. Therefore, such 
a PG should be detectable. On the other hand, the long 
dashed curve corresponds to a PG with a t s at the right 
end of the relevant time interval. The Ly-a luminosity of 
such a PG has just reached detectable values and will rise 
further. 

The maximal number of PGs in the volume A V c which 
we can expect to be bright enough in Ly-a is then given 
by the integral over the mass function of PGs with lower 
integration limit M m i n multiplied by AV C : 



$(M, z )dM AV C 



(3) 



M„ 



where the mass function <I>(M, z )dM gives the comov- 
ing density of objects with mass M £ [M, M + dM] at the 
observed redshift zq. 

But not all of these galaxies are in the PG phase at 
the observed epoch to- Some of them already left the PG 
phase or others enter this phase at a later epoch, so that 
they are dormant at to- The situation is illustrated in Fig. 



^ The Ly-a bright phase during which a PG will be de- 
tectable in Ly-a is constrained to a certain time interval 
[t s + Ati(M),t s + At 2 {M)\ which is defined by 



for all te [t s + At x (M) , t s + At 2 (M)] 



(4) 



The PG is observable in Ly-a at zq — z(to) only if to G 
[t s + Ati(M), t s + At2(M)]. This translates to a constraint 
for t s : 

t s € [t - At 2 (M),i - Ati(M)] (5) 

Only PGs with M > M mm which fulfill © will be de- 
tectable in Ly-a at a redshift of z (Fig.QJ. 

The next building block of our model is the 'history 
of galaxy formation', that is the distribution of ignition 
times t s which we assume to be coupled to the formation 
time of the halo. Galaxies show a substantial age spread. 
This translates to different t s for different galaxies. Finally 
our model should end up with the number of spheroids 
we observe today, e.g. with the current mass function of 
spheroids <&(M, 0). To take into account the 'history of 
galaxy formation', we introduce a distribution function 
PM(t)dt, which gives the fraction of current spheroids with 
mass M which formed and started their first star forma- 
tion during the time interval [t,t + dt\. The index AI in- 
dicates that Pjvf(i) will in general depend on the mass of 
the objects. We normalize Pm(^) according to 

t(z=0) 

P M (t)dt = 1 (6) 



t(z = 0) is the present age of the universe. 
$(M, 0) Pm (t)dtdM gives the comoving number density 
of spheroids with mass in the mass interval [M, M + dM } 
collapsing and starting their first star formation during 
the time interval [t, t + dt]. The number of PGs per solid 
angle Auj with a detectable Ly — a flux above Su m is then 
given by 



N, 



Ly a >L„ 



(7) 



lM mM {Ct^M) HM, 0)P M (t)dt} dM AV c (z ) 
Some remarks on equation (7): 

(i) In this form, expression (7) is quite general but note 
that the integration limits to — At 2 (Af) and to — 
Ati(M) are determined by the duration of the phase 
in which the PGs of a certain mass are bright enough 
in Ly-a to be detectable. They depend strongly on the 
shape of the function f(t — t s ). 

(ii) The mass function (for spheroids) at redshift z(t) is 
given by 



tO) 

$(M,z(t)) = $(M,0) / P M {t)dt 



(8) 



With the normalization © we ensure that the number 
of spheroids which form in our model is the number we 
observe today. 
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(iii) Note that with the normalization 10 we assume im- 
plicitly that the spheroids which form at high red- 
shift remain until today and therefore we neglect that 
some spheroids might merge after their formation. This 
means that our approach might underestimate the 
abundance of PGs at intermediate z. This could in 
principle be dealt with by adjusting the right hand 
side of © to values > 1. Because the merger rate 
might depend on the mass of the objects, these val- 
ues would clearly be mass dependent. However, in this 
paper we keep the normalization constraint © for sim- 
plicity and therefore neglect merging of spheroids after 
their formation. 

4. Discussion of two special cases 

In order to make expression (JJJ more transparent, we 
choose two special functions for Pu(t) and L^f^(t — t s ). 
First, we assume for Pjvr(t) a delta- function Pu{t) = 
8(t — ti n ). That is, all spheroids have the same formation 
and ignition epoch t s — ti n and shine in Ly-a simulta- 
neously (independent of their mass). If we assume that a 
fraction e of the stars of a spheroid are born in the first 
starburst with a constant SFR over a period At (that is 
SFR = eM},/ At, where Mf, is the stellar mass w baryonic 
mass of the spheroid today) and that the Ly-a luminosity 
is proportional to the SFR, we get for the Ly-a luminosity 



W*^£r if te [*»,*» + At] 

otherwise 

-(Mi 



(9) 



With this choice for P M (t) and L^(t - t s ) Q gives 

Ndet(z ) = 



otherwise 



(10) 



This is the formula Baron & White l|1987fl used for their 
calculations. 

They took t in = \t coU and At = ^t coU , where t coU 
is the duration of the collapse associated with a uni- 
form spherical perturbation of the same initial mean den- 
sity as the protogalaxy . They varied z co u between 6 
and 1.5 and took qo = 0.5 and qo = 0.05. This cor- 
responds to a variation of At between 0.6 and 6 Gyrs. 
For the luminosity function &(L) of present-day galax- 
ies they took a Schechter function with the parameters 
a = -1.25, L* = 1.6 x 10 10 /iq 2 L and $* = 1.2 x 10~ 2 ^ 
Mpc~ 3 . Assuming (as we do) that one only can observe 
the ancestors of present-day ellipticals and the bulges of 
present-day spirals during their violent formation pro 
cess, they reduced $* by a factor of 3 to 0.4 x 10~ 2 
Mpc -3 . They used a mass-to-light ratio of 6.6 to con- 
vert the present day luminosity function to a mass func- 
tion <!>(M, zq). To convert star formation rates into a Ly- 
a flux, they used for the constant kj^ ya in JjJJ the value 

1 Note that the collapse of a spherical perturbation starts at 
the redshift z coU = z(t con ). 



h 3 
"o 




o — 
1c 



18.5 -19 -19.5 -20 

log(S ,. m [W/m 2 ]) 



Fig. 2. Most optimistic (solid line) and pessimistic 
( dashed line) number density predictions for Ly-a emitting 
PGs for zq = 4.8 and Az = 0.1 according to the parame- 
ters of Baron & White (1987) together with the early upper 
limits of the Palomar Fabry-Perot survey of Thompson et 
al. a]) (triangles) and more recent searches by Maier 

et al. ( small arrows ) and Rhoads et al. ( open 
star with arrow). 



'zbw ■= 0.25 x 10 35 Sto-^ — r . Note that this value is by 
m- Moyr- 1 J 

a factor 4 smaller than the value deduced from local galax- 
ies with the assumption of a standard IMF and CASE B 
recombination (Kennicutt I1983fl . In this respect, Baron 
and White were conservative. Fig. [21 shows the results we 
got from l|10|) for zo = 4.8 and Az = 0.1 with the parame- 
ters of Baron & White together with the early limits from 
the Palomar Fabry-Perot survey for Ly-a emitting PGs 
(Thompson et al. 1995a, 1995b) and more recent limits 
from CADIS (Maier etal. fMty and LALA (Rhoads etal. 
2003). Even the least optimistic predictions by Baron & 
White are in obvious conflict with the results of the latest 
surveys. 

Alternatively, let us assume that the galaxies do not 
form and start at the same time with their Ly-a bright 
PG phase. A simple way to approximate this is to assume 
that the formation and ignition times t s of the galaxies 
are distributed equally over a certain time interval. So the 
delta function for Pm (t) has to be replaced by the function 



P M {t) 







otherwise 



(ii) 



Furthermore, we assume that the genuine Ly-a bright PG 
phase only lasts for a limited period, e.g. Athya = 0.2 Gyr 
as discussed above: 



kr €Mb 
^Lya At 



for te[t s , 
otherwise 



t s + At 



Lya\ 



(12) 
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Table 1. Factors AtL ya / Atp for different Zi n and 
and At Lya = 2 x 10 V (for H = 70km/s Mpc" 1 ). 



%out 





Zout 


f2_M 




Atp 


20 


5 


1.0 


0.0 


0.37 






0.3 


0.0 


0.23 






0.1 


0.0 


0.16 






0.3 


0.7 


0.20 






0.1 


0.9 


0.12 


20 


2 


1.0 
0.3 
0.1 
0.3 
0.1 


0.0 
0.0 
0.0 
0.7 
0.9 


0.12 
0.08 
0.06 
0.07 
0.04 



With HTJ and (JEJl we get 

N L y a>Lmzn (ZQ) = 

with 

St < At Ly - a 



(13) 



(14) 



Thus, in comparison with (|10f> the numbers are at least 
reduced by the factor AtL y _ a / Atp. In Tabled we show 
this factor for different z in , z out (z in and z out are the red- 
shifts, corresponding to t in and t out in JTT}) and a AtLya 
of 2 x 10 8 yr for different values of the matter density Q.m 
and the vacuum energy density fly- Obviously the fac- 
tors by which the numbers are reduced can be quite large. 
This simple example demonstrates that the predictions 
e.g. from Baron & White l|1987f) are likely to be much too 
optimistic and could well over-predict the abundance of 
Ly-a bright galaxies by a factor of 10 or so. 



5. The detailed model 

The above examples use oversimplified assumptions. A 
proper model instead should aim for more realistic func- 
tions L L A p a (t), P M (t) and $(M,0). We deliberately follow 
a phenomenological approach: We try to parameterize the 
situation with reasonable functions and try to fix the pa- 
rameters using both the present day galaxy population 
and some results from high redshift observations. Note 
the close similarity to the "semianalytic" models of galaxy 
evolution (see e.g. Somerville etal. |2001). 



5.1. The 'history' of galaxy formation 

In the previous section we approximated the function 
Pm(£) by assuming that the ignition times t s of the 
galaxies are either fixed at a certain time or equally dis- 
tributed over a time interval. Here we will estimate the 



function Pjvf(i) from the paradigm that galaxies are as- 
sumed to arise from peaks in the density field S(x,t) = 
{p{x,t) — p)/p. In the following we will omit the spatial 
coordinate x. 

Consider the density field S(M, ti) at an initial time 
ti smoothed with a box containing the mass M and let 
cro(M,ti) be the rms variation of this smoothed density 
field. We define the dimensionless fluctuation amplitude 



u{M,U) 



S(M,U) 
a (M, ti) 



(15) 



In the linear regime the fluctuation amplitude grows pro- 
portionally to the linear growth factor D(t). The density 
field at an epoch t > t iy therefore, is given by 



S{M,i) = S(M,U)D(t)/D(ti) 



(16) 



The density contrast grows according to this simple 
relation until non-linear effects become important and the 
region ceases to expand, turns around and collapses to 
form a virialized halo. According to the linear relation l|16fl 
the density contrast would have reached a critical value 
S c at the time the halos form. S c can be estimated from 
the evolution of an isolated spherical over-dense region 
(see e.g. Padmanabhan (1993)). The value is of order one 
and depends weakly on Qm- The exact value will not be 
important for our considerations. According to this rule, 
peaks in the initial density field (at epoch ti ) with a peak 
hight of v(M,U) = 5 c D(ti)/(a (M,ti)D(t )) will collapse 
at the epoch to- We define the function 



Ui{M,t) 



S C D(U) 



a (M,ti)D(t) 



(17) 



We are interested in the number density dNe orjx Jt) of 
peaks which collapse in a certain time interval [to, to + dt\. 
It depends on the distribution of the peak hights, p{v). 
p(y)dv gives the probability that a peak has a v- value in 
the interval [v, v + dv\. p{v) translates into a distribution 
of collapse times. dN^ orm (t) is given by 



dN, 



form 



p{ui{M,t)) 



dvi 



dt 



dt 



(18) 



from which we obtain of the distribution in collapse times 
dvi 



P M {t) =p(ui(M,t)) 



dt 



(19) 



The distribution of collapse redshifts is then given by 



P M (z)=p(H(M,t)) 



dvi 



dt 



dt 

dz 



(20) 



If we assume that we can approximate the fluctuation 
spectrum by a power-law with index n, we can write the 
function Vi (M, t) 



MM,t)=[ — 



M 
M* 



i + 3 
6 



D{U) 
D(t) 



(21) 
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Fig. 3. The distribution function p(v) for 7 = 0.58 cor- 
responding to an effective power-law index of the power 
spectrum of u = —2. The solid curve shows the result 
of the exact calculation and the dashed one is calculated 
with the fitting formula given by BBKS. 

M* is the characteristic mass scale, which collapses at the 
time t*. We treat for mass M* as a free parameter of 
our model. In this way we incorporate uncertainties in S c 
and the absolute normalization of the fluctuation power 
spectrum. 

If the density fluctuation field is Gaussian, we can cal- 
culate the peak distribution p(v) with the formalism given 
by Bardeen et al. I|19861 BBKS). p(v) depends on the 
width of the power spectrum under consideration. This is 
described by the BBKS parameter 7 (7 = 1 corresponds 
to power at a single wavelength only; lower values indi- 
cate a larger range of wavelength). In the special case of 
a power-law spectrum with a Gaussian filtering 7 is given 
by (BBKS) 



Assuming a cold dark matter (CDM) model with adiabatic 
initial density fluctuations with a power spectrum ~ fc, 
the power spectrum evolves as the universe expands into 
a power spectrum with effective power-law indices ranging 
from -3 on small scales up to 1 on large scales (see BBKS). 
On galactic scales the index is n=-2, which is the value we 
assume in the following. With Gaussian filtering l|22|l gives 
7 = 0.58. 

Since, in general, p(v) cannot be written in closed form, 
BBKS give a fitting formula which is accurate enough for 
our purpose. Fig. |3] shows p(v) for 7 = 0.58 calculated 
exactly and with the fitting formula provided by BBKS. 
We normalize Pm(^) according to equation ©. Fig. 0] 
shows the distribution Pm(0 of the collapse times for 
three different masses 0.1 * M*,M* and 5 x M*. Fig. 
[5] shows the corresponding Pm(z) curves, t* was cho- 
sen in such a way that PM„{t{z)) has its maximum at 
zmax = ^(tniax) = 6. The curves illustrate again the cru- 
cial point of 'bottom up' hierarchical structure formation: 



t [Gyr] 

Fig. 4. Distribution of collapse times for the three differ- 
ent masses 0.1 x M» (dotted-dashed line), Af» (solid line) 
and 5 x Af» (dashed line). 

PM(t(z)) with M < M* peaks at higher redshifts than 
Pm, (t(z)) and PM(t(z)) with M > M* peaks at lower red- 
shifts than PM,(t(z)). Furthermore, the scatter (width) of 
the collapse times is smaller for small mass objects which 
form at high z than for massive objects which form later. 
Note that the peak formalism is not a strong theoretical 
motivation for the shape of our function Pm(^)- By tak- 
ing z m ax as a free parameter of the model (see above) we 
are not strictly applying the peak formalism which would 
give a definite result for the mass function of halos at a 
certain redshift once the power spectrum is known. Here 
we only use the peak formalism to get an estimate for the 
form of Pm(^s)- Furthermore remember that t s is not the 
formation time of the halo but the time when the first 
star formation starts whereas in the mass function which 
would be predicted directly by the peak formalism the 
variable would be the time of collapse t co i of the halos. 
By using z max as a free parameter we leave open how ex- 
actly the collapse time and the ignition of star formation 
in halos are correlated. Comparison of our model with ob- 
servations may in the end determine wether our choice for 
PAi(t s ) according to the paradigm of hirarchical structure 
formation was a good choice or not. However, because of 
the simplicity of our model, we could easily modify the 
function Pm(^s) to take into account new results such 
as that bigger galaxies may have earlier star formation 
epochs as suggested by Heavens et al. • 20(51$. This may 
be the subject of a following paper. 

5.2. Development of the Ly-a luminosity with time 

The assumption of a constant star formation rate and Ly- 
a luminosity of a PG as in l|12H is a crude simplification. 
Detailed numerical simulations of the formation of galax- 
ies (e.g. Steinmetz Q^S3) show that the SFR of PGs start 
with very low values, increase rapidly and reach a peak 
after a few 10 8 years. The metallicity of the early bulge 
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Fig. 5. The same curves as in Fig. but with the time 
axis transformed into redshift. 

increases rapidly in this phase and soon reaches values 
around 1/10 of the solar value. According to our definition 
this limits the genuine Ly-a bright PG phase to a short 
time period At^ya of the order 10 8-9 years. After this time 
the metalicity reaches values which inevitably lead to sig- 
nificant dust formation combined with an absorption of 
the Ly-a flux. Although the escape of Ly-a photons from 
star forming galaxies is a rather complicated problem de- 
pending on several conditions like the composition of the 
interstellar medium (e.g. multiphases, see Neufeld ©91), 
or gas flows in the interstellar medium, see Kunth et al. 
(1999) , we come back to this points in section 7), we 
simplify and assume that after a certain time period the 
Ly-a emission is no longer proportional to the SFR and 
decreases mainly due to the increasing dust content. We 
approximate this behavior by a Gaussian for f(t — t s ) (see 
PI ) with a width a^ya which remains a free parameter 
of our model: 



(t — t S —2<TLy 

72 



2a 







-| ioxt>t s 
otherwise 



(23) 



Later on, we will see that the observed number densities of 
Ly-a emitting PGs constrains (JLya- The typical length of 
the Ly-a bright phase could be described by the FWHM 
AtLya = 2crL yQ V21n2 of J2H1). As already expressed in Q), 
we assume that the star formation rate and therefore the 
Ly-a flux during this first phase of star formation is pro- 
portional to the baryonic mass Mb content of the galaxy. 
We further assume that Mb is proportional to the total 
mass M of the object according to M& = j^M (f^ is the 
density parameter of the baryonic mass in the universe). 
Furthermore, we assume that the star formation rate and 
accordingly the Ly-a luminosity scales proportionally to 
l/t s . This is motivated by the spherical collapse model 
(see e.g. Padmanabhan . 1 1 993[) . A perturbation which de- 
couples from the cosmic expansion at the epoch tmax will 



collapse (to a black hole) at the epoch t s = 2 x t max - The 
duration of the collapse process is therefore St — 0.5 x t s . 
Hence, objects which collapse earlier collapse faster and 
we assume that the star formation rate is higher as well. 
Since we do not attempt to incorporate the complicated 
processes leading to primeval star formation and a certain 
IMF into our model, we summarize this in the free pa- 
rameter e, which describes the overall scaling of the star 
formation with respect to Mb. As in our simple model in 
section 4 (see equation ©) , e can be interpreted as the 
factor of proportionality which determines the fraction of 
the stars of the spheroid which are born in the first star- 
burst, although this is not strictly true as the Gauss func- 
tion l|23(l describes the combined effect of star formation 
and dust production on the Ly-a emission. If all stars are 
produced in a very short star burst, e could be well above 
1. In this paper we fixed e to one. This means that the 
Ly-a emission of a galaxy at its maximum corresponds to 
a SFR of SFR = 

5.3. Mass function 

We assume that 

— Ly-a emitting PGs are the precursors of todays 
spheroids 

— the baryonic mass we see today in these spheroids is 
proportional to the baryonic or total mass of the cor- 
responding PG at the time of their first star formation 

— subsequent merging processes or gas accretion does not 
destroy the luminous part of the average galaxy, which 
forms todays spheroid. However, merging of the PG 
into a larger halo (e.g. a cluster) may strip gas and/or 
the original dark halo of the galaxy and gas accre- 
tion might lead to the formation of a disk. Additional 
star bursts might add further stars to the spheroidal 
system, but the added baryonic mass is assumed to 
be proportional to the baryonic/total mass initially 
present at the time of the first star burst. 

In equation 123fl we therefore need the todays (bary- 
onic) mass function of spheroids ^spheroids 0^&> 0)- We 
derive the local bulge mass function from the type depen- 
dent luminosity functions determined from the CFA red- 
shift survey by Marzke etal. (1994) together with values 
for the type-dependent ratio of the bulge to disk lumi- 
nosity from Simien & De Vaucouleurs |1986|l and a con- 
stant (baryonic) mass to light ratio for spheroids set to 
1OM /L O . 

Marzke et al. I|1994f) described the luminosity functions 
for all Hubble types T by Schechter functions 



§ T {L)dL = [ — j cxp 



L 



L 



d[ — ) (24) 

J T / V . 

with the parameters listed in table 

We write the total luminosity L to t of a galaxy as the 
sum of the bulge luminosity Lb and the disk luminosity 
Ls- 



Ltot = Lb + Ls 



(25) 
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Table 2. Schechter function parameter for the different 
Hubble types according to Marzke et.al. (1994) 



Hubble Typ T 


a 


L* 


$* 






[W 9 h- 2 L Q ] 


[10- 3 /igMp C - 3 ] 


E T G [-7, -4J 


-0.85 


7.66 


1.5 


SO T6[-3,0] 


-0.94 


4.88 


7.6 


Sa-Sb T£[l,4] 


-0.58 


4.79 


8.7 


Sc-Sd T G [5, 7] 


-0.96 


5.2 


4.4 


Sm-Im T > 7 


-1.87 


5.1 


0.6 



Table 3. Effective bulge mass-to-light ratio Tjf for the 
different Hubble types. 



Hubble Typ T 


[M Q /L e ] 


e re [-7, -4] 


10 


SO T G [-3, 0] 


5.4 


Sa-Sb Te [1,4] 


2.8 


Sc-Sd T G [5, 7] 


0.48 



Furthermore, we denote with the ratio of the bulge 
luminosity to the disk luminosity: 

IT ■= ^ (26) 

The bulge mass of a galaxy of Hubble type T is given by 

M B = T B 1 =: T e J f L tot (27) 

1 + 7 T 

Tb denotes the mass-to-light ratio of the bulge compo- 
nents of present day galaxies, which is relatively inde- 
pendent of the galaxy type. We took = 10 M /L Q . 
Eq. I|27|l defines T^/^, a type dependent "effective mass- 
to-light ratio". Table |21 lists the values Tyf^ for different 
Hubble types calculated with the type dependent jt val- 
ues of Simien and DeVaucouleurs (1986). 

M B (T) = T e J f L* T (28) 

With (|28|l the luminosity functions l|24|l transform into the 
type dependent bulge massfunctions 

$ T (M B )dM B = 

(29) 

^(sfjrj) ex p{-A7ffr)} d (l^fr)) 

The total (baryonic) spheroid mass function of present day 
galaxies is given by the sum over the mass functions of the 
different Hubble types 

$tot (M B ) dM B =J2® T ( Ms ) dMe ( 3 °) 

T 

We will use as an approximation for the mass function 
$(Mf,,0) in equation J2J. 



5.4. Summary of the model 

We constructed a phenomenological model to calculate the 
expected surface density of Ly-a emitting, high redshift 
young galaxies. According to Q we have to take into ac- 
count the evolution of the Ly-a luminosity as a function 
of time and mass of the galaxy and the galaxy formation 
history. We estimated these functions with the following 
assumptions: 

— We assume that Ly-a emitting PGs are the precursors 
of the present day spheroids (=bulges & ellipticals) in 
the phase of their first star burst. 

— The first star formation starts in a dust-free environ- 
ment. The enrichment of the 1GM with dust in the 
first hundred million years of star formation leads to a 
rapid attenuation of the Ly-a emission. The time evo- 
lution of the Ly-a emission is described by a Gaussian 
with FWHM At^ya, which gives the typical duration 
of the Ly-a bright phase. 

— The SFR and accordingly the Ly-a emission is as- 
sumed to be proportional to the baryonic mass Mf, 
of the PG. Furthermore, the SFR is assumed to scale 
with j-, where t s is the collapse time: Comparing ob- 
jects with the same baryonic mass, the ones which form 
their stars at higher redshifts have higher star forma- 
tion rates. 

— We get a first guess for the functional form of the dis- 
tribution of 'ignition times' Pm(^s) by the distribution 
of peak hights p[y). However, this is a mere ad-hoc as- 
sumption with no strong theoretical motivation. We fix 
the parameters M* and t* so that P]\i(t s ) peaks at a 
certain redshift z ma x for M = M* = 4 x 10 10 M o , 
which corresponds roughly to the (baryonic) bulge 
mass of our Milky Way. Zmax is a free parameter of 
our models. 

— We normalized Puifs) so that our model ends up at 
z = with the mass function of spheroids we see today. 
This assumes that the spheroids which formed in the 
early universe stay as entities until today and that their 
current baryonic mass is proportional to the baryonic 
mass at the time of their first star formation. This does 
not exclude later merging processes of the dark matter 
halos and e.g. accretion of gas, which may lead to the 
formation of a disk. 

An overview of parameters of our model is given in Table 
0J Note that we keep those parameters fixed which are 
determined by observations of the local universe or theo- 
retical arguments. The free parameters of our model have 
to be determined by the observational statistics of Ly-a 
bright PGs. This will be attempted in paper II. Here we 
will only point to the most obvious constraints (see section 
7.2) and use them to define a "basic model" . In Table 4 the 
parameter values of this "basic model" are summarized. 

6. Predicted abundance of Ly-a galaxies 

Fig.H3shows the expected number density of Ly-a emitting 
PGs as a function of zmax- The different panels show the 
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Table 4. Parameters of the model and their value range explored in this paper. The values of the "basic model " are 
chosen in such a way that the model is in agreement with recent observations (see section 7.2 ) 



Parameter 


Value range 


Explanation 


"Basic Model" 




n M = 1.0, 0\ = 0.0 

n M = 0.3, n A = o.o 
n M = 0.3, n A = 0.7 


flat CDM (fCDM) 
open CDM (oCDM) 
A CDM 


n M = 0.3, O a = 0.7 


h 


0.7 


Hubble paramter Hq = /ilOOkm /s Mpc~ x 


0.7 


At Lya = 2CT ij/a \/21n2 


0.1 - 2 Gyr 


duration of the Ly-a bright phase 


0.35 Gyr 


e 


e ee 1 


star formation effiency 


1 


n 


-2.0 


power law index of the power spectrum 


-2.0 


Zmax 


3-25 


redshift, at which the distribution of 
'ignition times ' Pjf(t,) peaks for 
M = M,=4x 10 10 M Q 


3.4 


a,L*, $*, T e " 


see Table 2 and Table 3 


parameters describing the luminosity function and 
the mass to light ratio for todays spheroids 


see Table 2 and 
Table 3 





i i i i i i i i i i i 
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Fig. 6. Expected number density per deg 2 and Az = 0.1 of Ly-a emitting PGs with a minimum detectable Ly-a flux 
of Sum = 3 x 10~ 20 Wm~ 2 . AtLya was fixed to 0.35 Gyr according to our "basic model". The different panels show the 
results for different observing redshifts zq. The three different curves in each panel correspond to different cosmologies: 
The solid line corresponds to Am = 0.3, Qa — 0.0, the dashed line corresponds to Qm = 1.0, f^A = 0.0 and the dashed 
dotted line corresponds to Om = 0.3, Q\ = 0.7. 
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results for different observing redshifts z . Note that the 
curves for the different redshifts look very similar. They 
rise steeply at small z rnax < z to a maximum and then 
fall off slowly at high redshifts z max > Zq. This behavior 
can be understood by Fig. [7| which shows Pm( z ) for ellip- 
tical galaxies with M = M* and various z max — 3, 6, 9 
and 12. The two vertical solid lines mark the interval (J3J) 
for zq = 3.5. The number of observable Ly-a emitting PGs 
is proportional to the shaded area under the Pm{z) curve 
in this interval. The observed number density is governed 
by two effects: For low z max the maximum of the Pm(z) 
curve lies on the left hand side of the interval (JSJ, that is 
at redshifts below zq (see the solid curve in Fig. UJ). With 
decreasing z max the width of the Pm(0 curve increases 
rapidly. This corresponds to an increasing spread of the ig- 
nition times, that is an increasing Atp in l|ll|l (see also the 
discussion there). Taking into account the normalization 
constraint JBJ), the overall amplitude of Pm(z) decreases 
rapidly with decreasing z m ax- This explains the steep de- 
crease of the curves in Fig. for decreasing z rnax < z$. 

If z max lies at higher redshift than the interval © 
we face the following situation: With increasing z max 
the overall amplitude of Pm(z) increases rapidly (see the 
dashed, dotted and dashed dotted lines in Fig.[7|).The rea- 
son for this is that the width of the curves Pm(^) decreases 
with increasing z max . However, because the total number 
of galaxies which form is conserved (see again ©) the 
total area under the Pm(^) curve does not change much. 
Thus, the peak amplitudes of Pm(^) and Pm(z) increase. 
This almost completely compensates for the effect that for 
Zmax > zq the observed redshift interval (see Fig. d) sam- 
ples lower and lower relative levels of the Pm(z) curve. 
If one does not consider one specific mass only but the 
full range of masses M > Mu m , the situation becomes 
somewhat more involved but the basic explanation for the 
behavior of the N(z max ) curves remains valid. 

The lines in the different panels of Fig. El reach roughly 
the same maximum value N max , almost independently of 
the observation redshift zq. For instance, for £l m ,Q\ = 
0.3,0.7, we find 1200 < N rnax < 1600/D°/(Az = 0.1) out 
to redshift zq = 9.3. This can be understood in the follow- 
ing way: (1) Because we assumed that the star formation 
rate and accordingly the Ly-a flux scales proportional to 
l/t s (see eq. (|23|l '). galaxies of a given mass are intrinsi- 
cally brighter in Ly-a if they form at higher redshift. (2) 
As discussed above, the function Pm(^) is narrower if it 
peaks at higher redshifts (see Fig. 01 and Fig. and thus 
exhibits a higher peak amplitude. 

(3) In addition one should note that the comoving vol- 
ume AV(z) is almost independent of z between z ~ 3 and 

(1) and (2) lead to an increase of PGs above a fixed 
mass limit with increasing z$. On the other hand, increas- 
ing zq means increasing the luminosity distance of the 
observed objects, and thus requires higher masses for PGs 
which are still observable above the detection limit Su m . 
Since both t(z) and Dl(z) directly depend on the geom- 
etry of the universe, the fact that N max is almost inde- 



K) 




2 4 6 



Fig. 7. Pm(£) for elliptical galaxies with M — M* and 
different z max — 3, 6, 9 and 12. The two vertical lines mark 
the interval (JSJ) for zq — 3.5. The dashed vertical line very 
close to the left border of the interval © marks the ob- 
serving redshfit zq — 3.5. 

pendent of zq is not a coincidence but rather a genuine 
property of our model. 

After to z max , the least known free parameter of our 
model is the duration of the Ly-a bright phase, AtL ya . 
Actually, it depends on the detailed astrophysical condi- 
tions during the formation of the first generation of mas- 
sive stars in PGs. Most notable are the initial mass func- 
tion (IMF) which determines the rate of heavy element 
production, the feedback processes which enrich the inter- 
stellar gas, the cooling of the hot gas phase, the topology 
of the star forming regions and many more. Therefore it 
is essential to understand how AtL ya = 2\/2 In 2 ■ <7L ya 
influences the predicted number of Ly-a galaxies in our 
model. 

The panels in Fig.|Slshow the expected number density 
of Ly-a emitting PGs as a function of ui, yol for a survey 
flux limit = 3 x 10~ 20 W/m 2 . From equation we 
get for the length At of the time interval © 



At = At 2 - Ah = 2a Lya J2ln(— : ) (31) 

First, we consider the case cfLya << t(zo), z max > zq. In 
this case the time interval (0 ends shortly before t(zo) 
and the maximum of the Pm (t) curve is located at higher 
redshifts than the redshift corresponding to (to — A^) (see 
Fig. 0). An increase of (Jhya has two effects which both 
increase the expected number density of Ly-a emitting 
PGs: 

— The interval 10 over which Pm(^) is integrated in- 
creases proportional to OL ya 

— The interval (0) moves to the right, that is nearer to 
the maximum of Pm(^)- 

If z max < zq and <7L ya << t(zo) the increase of the in- 
terval (5) proportional to a^ya gives an increase of the 
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expected number density, too. Accordingly, all curves in 
Fig. IHI have a positive slope at small aLya- On the other 
hand, if o- Lya is large, the interval JSJ moves to earlier 
times before Pm if) reaches its maximum and where Pm if) 
is a steep function of t. This leads to a decrease of the ex- 
pected number density of Ly-a emitting PGs when further 
increasing cri,y a . 

Finally, we show the expected abundance of Ly-a 
galaxies as function of the survey limit Su m and observed 
redshift zo (see Fig. EJ). For the purpose of this plot we 
fixed the duration of the Ly-a bright phase to the value 
of our basic model: AtL ya = 0.35 Gyrs (see section 7.2). 
To illustrate the dependence on z rnax we display the cu- 
mulative number density of Ly-a emitting PGs for three 
different values of z max = 3.4, 6.0 and 10.0. It is obvious 
from Fig. [5] that for z max < 10 a set of optical surveys 
for Lyman-a (zq = 3.5 corresponds to A(Lya) = 547 nm, 
zq = 6.6 to 924 nm) could be sufficient to determine z max . 
However, when z m ax is located beyond z max = 10, only 
the inclusion of a deep near infrared emission line survey 
(e.g. in the J band aiming for Ly-a at zq ~ 9.3) would 
be conclusive. Such a survey seems not feasible from the 
ground and would have to wait for the Next Generation 
Space Telescope. 

7. Discussion 

After we outlined the generic properties and predictions of 
our model, we discuss several aspects which relate to cur- 
rent and future surveys for Ly-a emitting primeval galax- 
ies. 

7.1. Optimum survey strategy 

First, we consider what conclusions about an optimum 
survey strategy for Ly-a galaxies can be drawn from the 
predicted number densities (Fig.[5|): At bright detection 
limits (e.g. log(Sii m ) ^ —19.5 for z < 6 ) the curves 
are steep due to the exponential fall-off of the underlying 
luminosity function (i.e. underlying mass function (j20l) 
towards high luminosities. Here it will be more useful to 
improve the detection limit of a survey than to enlarge 
its area. The opposite is true at very faint detection lim- 
its (e.g. log (Sii m ) ^ —20.5 for z < 6). A survey which 
reaches such a depth will benefit more from an increase in 
area than from pushing the limits deeper. 

Formally, one might quantify the merit of a Ly-a sur- 
vey by the total number Npq of Ly-a emitting PGs which 
could be found by spending a given observing time t b s at 
a given telescope. 

Using the observing time t obs to increase the observed 
area Au on the sky gives 

Ap G oc Aw oc t obs 

=> log(Ap G ) = log(t o6s ) + const. (32) 

aiog(jVpG) _ i 

^ d\og(t obs ) - 1 
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Fig. 10. Comparison of our "basic model" with observa- 
tional data. Big open star: LALA survey of Rhoads etal. 
(2003); filled square: Hu etal. (1999); small stars and ar- 
rows: CADIS Maier etal. I|2003fl : open boxes: Hu etal. 
(1998);open circles: Kudritzki etal. (|2000[> : filled star with 
arrow: LALA survey of Rhoads etal. U003 for zq = 4.8; 
solid line: zq = 3.5; dashed line: zq = 5.7; dotted dashed 
line: z — 4.8; 

On the other hand, using observing time to improve the 
detection limit <Sii m : 



Slim oc 



(33) 



=> Iog(Su m ) = -^og(t obs ) + const. 
If we take into account that log(Afpc) is a function of Sii m 

log(A^ PG ) = /(log(5 lim )) = f(- 1 -lo g (t obs )+ canst) (34) 



we find: 

dlog(JVpG) 

d\og(t obs ) 



-i/'(log(S lim )) 



(35) 



Notice that /' < everywhere (cf. Fig.^J). Comparison 
of (|32|l and Ij35(l shows how the optimum survey strategy 
should be chosen: If /'(log(Sii m )) < —2 one should use 
additional observing time to increase the integration time 
in order to reach fainter flux limits. If /'(log(5n m )) > —2 
additional observing time should be used to increase the 
survey area Alu on the sky in order to increase the number 
density of detectable PGs. 

7.2. Comparison with observed abundance 

We have chosen the parameters of our "basic model" in 
such a way that our model fits the observed surface density 
of Ly-a emitting PGs at redshifts zq = 3.5 and zq = 5.7. 
Fig. E3 shows the result. If we keep e fixed to 1 the main 



Thommes & Meisenheimer: Primeval Lyman-a galaxies 



13 




i 1 1 1 1 i 1 1 1 1 i 1 1 1 1 - 

4. 







-1 -0.5 0.5 



4 



z =9.i 



-0.5 
'Ly-a 



0.5 



log( (J [Gyr] ) 



-0.5 
log( CT Ly _„ [Gyr] ) 

Fig. 8. Expected number density per deg 2 and Az = 0.1 of Ly-a emitting PGs as a function of a^ya fsee !23[) for fixed 
Sum = 3x 10~ 20 W/m 2 . The cosmology is fixed to Qm — 0.3, Cl\ = 0.7. The solid curves correspond to z max = 3.4, 
the dashed curves to z max — 6 and the dotted dashed curves to z max = 10. The arrow marks the value ULy-a = 0.15 
of our "basic model" . 



free paramters of the model are z max and Ai£ y _ Q . These 
parameters seem to be strongly constrained by the obser- 
vational data at different redshifts. To produce Fig. 1 101 we 
varied the two parameters z max and Ai^y-a by hand in 
order to fit the data. It is non trivial to find a combina- 
tion of values consistent with all available data. However, 
as explained below, part of the difficulty may arise from a 
second Ly-a bright phase of PGs. 

FigEl shows the surface density of Ly-a emitters as a 
function of redshifts zq for different z max and a fixed de- 
tection limit. Fig. ll2l shows the same diagram for different 
detection limits log(£ Km [W/m 2 ]) =-21, -20, -19 but fixed 
Zmax = 3.4. At faint detection limits log(Szi m [W/m 2 ]) < 
—20.0 the number density peaks at redshifts zq > 3 and 
falls off slowly with increasing redshift. For z max =4.5 
and 5.5 the surface density of Ly-a emitters at a detec- 
tion limit of log(5Yi m [W/m 2 ]) < —20.0 is expected to be 
more or less constant in the redshift range zq = 3... 5 with 
values > 10 3 1/deg 2 / Az — 0.1, reaching its maximum 
at z w 4. For z max = 3.4 the surface density falls faster 
with increasing redshift zo > 3, especially at lower de- 
tection limits. However, at faint enough detection limits, 
log(Sz im [W/m 2 ]) w-21 (see Fig. U2J), the number density 
only changes moderately in the redshift range between 3 
and 5. At faint detection limits log(Szi TO [W/m 2 ]) < —20.0 
the expected number density of Ly-a emitting PGs is ex- 



pected to be significantly higher (> several 100 / deg 2 
/ Az = 0.1) out to very high redshifts of zq ps 8. This 
may open interesting prospects for the examination of the 
epoch of reionisation (see Haiman, 2002]and Cen, |2003). 
On the other hand according to Fig.^|the surface density 
of Ly-a emitting PGs is predicted to be a steep function of 
observing redshfit zq and detection limit at high detection 
limits log(S'/i m W/m 2 ]) > —19 which explains the failure 
of early surveys (see Pritchet et al. Q.994). 

As we determined the model parameters only with the 
data at zq = 3.5 and zq = 5.7 we get an independent pre- 
diction for the surface density of Ly-a emitting PGs at e.g. 
z = 4.8 (see Fig. llOf) . According to this prediction, the sur- 
face density at z = 4.8 e.g. at log(Su m [W/m 2 ]) = —19.5 
should be roughly a factor of 4 higher than at z = 5.7. But 
very recent data by the LALA survey (Rhoads et al. 2003 ) , 
CADIS survey (Maier et al. I2003|l and the Subaru Deep 
Field survey (Shimasaku et al. [2003 ) indicate that the sur- 
face density of Ly-a emitting PGs does not differ much be- 
tween z = 4.8 and z = 5.7 at log(S' ;im [W/m 2 ]) = -19.5. 
As we will discuss in detail in paper II, this can be under- 
stood in the framework of our models if z max is close to 
5 or 6 instead of 3.4 (cf. Fig.^2f° r the case z max — 5.5). 
However, the low absolute numbers observed at zq = 5.7 
and zq — 4.8 can then only be explained by reducing ei- 
ther e or o Lya considerably. As this will also reduce the 
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Fig. 9. Expected number density per deg 2 and Az = 0.1 of Ly-a emitting PGs as a function of the minimum detectable 
Ly Q flux of Sum according to our "basic model" parameters AtL y - a — 0.35 Gyr and Q,m = 0.3, Qy = 0.7. The solid 



curves correspond to z max = 3.4, the dashed curves to z v 



6 and the dashed dotted curves to z„ 



10. 



absolute numbers at zq = 3.5 a conflict with the obser- 
vations is unavoidable. This might be a hint that a large 
fraction of the Ly-a emitters at zq < 4 do not belong to 
the class of "primeval" Ly-a galaxies considered here, but 
are galaxies in a later state of evolution (see section 7.4) 

7.3. Comparison with other models 

Besides the calculations of Haiman and Spaans (1999| our 
calculations are the only ones that try to model the num- 
ber density of Ly-a emitters at high redshifts taking into 
account the formation of dust during the early phases 
of star formation and the history of galaxy formation. 
Haiman and Spaans deduced the mass function and for- 
mation history of haloes directly from the power spectrum 
with the Press-Schechter formalism, whereas we extrapo- 
lated the local luminostiy function of galaxies and their 
stellar content back into the past and only used the power 
spectrum and peak formalism to deduce a realistic distri- 
bution of formation times. Furthermore, while we use a 
phcnomenological approach to describe the modulation of 



the Ly-a emission by dust formation in the early phase of 
star formation, Haiman and Spaans used detailed Monte 
Carlo simulations of individual galaxies with a range of 
masses for the ionizing stars, dust content and inhomo- 
geneity together with the solutions of the radiative trans- 
fer problem for the Ly-a line in an inhomogeneous multi- 
phase medium (Neufeld, D-991). It is interesting that our 
model and the model of Haiman and Spaans both agree 
in predicting high surface densities of Ly-a emitters out 
to redshifts » 8. Furthermore, both models also predict 
that the surface density of Ly-a emitters as a function of 
redshift at faint detection limits (see Fig. I12|) should be 
rather flat in the redshift interval between 4 and 6. 

7.4. Lyman-break galaxies and second generation Ly-a 
emission. 

In our models we assumed that strong Ly-a emitting 
PGs are young spheroids during their very first phase 
of star formation in which dust does not play a crucial 
role. Subsequently the interstellar medium will be enriched 
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Fig. 11. Expected number density per deg 2 and Az = 0.1 
of Ly — a emitting PGs as a function of the observa- 
tion redshift zq for AtL ya — ^Ly-a 

= 0.35 Gyr 

(according to our "basic model") at a detection limit of 
^(■S'zimlW/m 2 ]) = —19.5. The three curves correspond to 
Zmax = 3.4 (solid curve) ,4.5 (dashed curve), 5.5 (dashed 
dotted curve). The cosmology is fixed to £1m = 0.3, 
ftv = 0.7 (according to our "basic model"). The sym- 
bols mark observational data: open box: Hu etal. ([1998) 
for log(Sn m [W/m 2 ]) = —19.5; open circle: Kudritzki etal. 
lEITO for log(S Hm [W/m 2 ]) = -19.5; small star: CADIS 
Maier etal. (|20Ti3|l for log(S Um [W /m 2 }) = -19.54. 
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Fig. 12. Expected number density per deg 2 and Az = 0.1 
of Ly — a emitting PGs as a function of the observation 
redshift zq for AtLya — 2oL y - a \/2\n2 = 0.35 Gyr and 
the detection limits of log{Su m [W/m 2 ]) = -21.0 (solid 
line), -20.0 (dashed line), -19.0 (dashed dotted line). z max 
is fixed to 3.4 and the cosmology is fixed to = 0.3, 
ily — 0.7 (according to our "basic model"). The symbols 
mark observational data: open box: Hu etal. (1993 for 
log(Su m [W/m 2 ]) = -19.1,-19,9; open circle: Kudritzki 
etal. lt2TW)ll for log(S Hm [ W /m 2 ] ) = -19.0,-19.8; small 
stars: CADIS Maier etal. lj2UU3|) for log(Si un [W /m 2 }) = 
-19.54,-19.32,-19.19. 



with metals very soon after the onset of star formation. 
Due to the ongoing dust formation, Ly-a photons will be 
more and more absorbed. The Ly-a flux of the PGs will 
decrease although the SFR may still increase and reach its 
maximum at a later time. A Ly-a dark phase follows, dur- 
ing which all Ly-a photons are destroyed due to resonant 
scattering in the dusty interstellar medium and during 
which the SFR reaches its maximum. 

At later times after the SFR reaches its maximum, 
model calculations (see e.g. Friaga and Terlevich (1999)) 
predict a phase in which strong outflowing winds build up. 
Once these winds have developed, a second Ly-a bright 
phase might develop. Complicated outflows of gas with 
high velocities v are indeed a common feature of LBGs 
(see Pettini et aL U998l and have also been found in nearby 
HII galaxies (Kunth etal. 1998). Large scale outflows not 
only explain that whenever LBGs show Ly-a in emission, 
this is shifted by up to « lOOOkm/s relative to the metal 
absorption lines but in addition explain their P-Cygni line 
profiles. 

The Ly-a line of objects in this second Ly-a bright phase 
should be shifted by the velocity of the outflowing wind 
relative to the metal absorption lines and should show 
a P-cygni profile. In principle high resolution high S/N 
spectra of the objects may allow us to disentangle Ly-a 
emitting objects which are in their first or second Ly-a 
bright phase. In our predictions for the number density 



of Ly-a emitters we took into account only the first Ly-a 
bright phase and neglected the possible second phase. At 
very faint detection limits and redshift z n < 4, Ly-a emit- 
ters from the second phase might contribute. However, 
because we fixed our parameters to the observed surface 
density of Ly-a emitters by Hu etal. (1998) and the new 
CADIS (Maier et al. 120031) results, which do not distin- 
guish between two Ly-a emitting phases of PGs, our value 
for AtLya might be representative of a (weighted) sum 
of the durations of both phases. In this, the duration of 
the second phase might have a lower weight, because this 
phase might be much fainter in Ly-a than the first one. 
Furthermore, one would expect that the relative "contam- 
ination" by Ly-a emitters in this second Ly-a bright phase 
increases with time. This might explain why Ly-a emit- 
ters at zo = 3.5 seem overabundant in comparison to a 
model which fits their density at redshifts zq > 4.5. 

8. Conclusions 

We presented a simple phenomenological quantitative 
model for the expected surface number density of high 
redshift Ly-a emitting galaxies. We assumed that ellipti- 
cal galaxies and bulges of spiral galaxies (which we call 
spheroids) formed early in the universe while disks were 
built up at a later stage. Thus, we identified the high red- 
shift Ly-a emitting PGs with these spheroids during their 
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first burst of star formation. One of the main assumptions 
of our model is that the Ly-a bright phase of this first 
starburst is confined to the first several hundred million 
years after the onset of star formation (duration: AtL ya ). 
We assumed an ad-hoc-function for the distribution of 'ig- 
nition times' with some motivation from the distribution 
of peak hights in the peak-formalism. In order to derive 
absolute number densities, we follow the method back- 
wards in time as pioneered by Baron & White (1987): The 
number of PGs that form in our model are normalized to 
the present (baryonic) mass function of spheroids. Using 
the surface density of Ly-a emitters detected by recent 
surveys at redshifts 3.5 and 5.7, we find that the Ly-a 
bright phase of primeval galaxies is very likely confined 
to a rather short period of < 0.5 Gyr after the onset of 
star formation. Our model predicts that the surface den- 
sity of Ly-a emitters with Ly-a fluxes <% m < 10 _20 W/m 2 
should be high (> several 100/D°/ (Az = 0.1)) out to very 
high redshifts of zq ~ 8. The substantial number of spec- 
troscopically confirmed high redshift Ly-a emitting ob- 
jects at redshifts zq > 5 (see e.g. Santos et al. 120041 or 
Malhotra and Rhoads 2004] and references therein) show 
that systematic searches for these objects are indeed suc- 
cessful. Now the main task for observers will be to quantify 
the selection effects and to separate the second generation 
Ly-a emitters. Together with our simple phenomenologi- 
cal model the observation of the Ly-a luminosity function 
at high redshifts (e.g. z — 5.7,6.6,9.3,12.6 as discussed 
above) may give interesting hints concerning the peak of 
galaxy formation activity (from z max ) and the duration 
of Ly-a bright phases of PGs (from AtL y - a )- As soon as 
this is accomplished our model will be able to pin down 
the formation of galaxies in term of the three parameters 
e, At^y^a and z max . It will be straightforward to test any 
physical model of galaxy formation with respect to this 
intuitive parameterisation of the observed abundance of 
PGs. 
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